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in 2+1- and 3- flavor QCD using 16 3 x 4 and 16 3 x 6 lattices and improved staggered p4 action. 
We also compare the renormalized Polyakov loop with the results of earlier studies. 
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1. Singlet free energies 

We did finite temperature lattice calculations in 3- and (2+1)- flavor QCD in the region of small 
quark masses. In the case of 2+1 flavor QCD we used two light quark masses m q = 0Am s and 0.2m s 
with m s being the strange quark mass. In the case of degenerate three flavors the quark masses were 
roughly 0.15/n $ and 0.3m s . Calculations have be done on 16 3 x 4 and 16 3 x 6 lattices. The lattice 
spacing and thus the temperature scale has been fixed using the Sommer scale = 0.469 fm [Q]. 
We used the interpolation Ansatz for the dependence of ro on the lattice gauge coupling j3 =6/g 2 
given in Ref . [[!]] . In our simulations we used the exact RHMC algorithm for the 2+1 flavor case 
while the standard R-algorithm was used in the 3 flavor calculations. Further details about our 
simulations can be found in Refs. ^ Qj. 

On the gauge configurations separated by 50 trajectories we have calculated the singlet free 
energy of a static quark anti-quark pair defined as 

exp(- J F 1 (r,r)/r+C) = -(TrW(r)W^(0)), (1.1) 
with W (r) being the temporal Wilson line. The above definition requires gauge fixing and we use 



the Coulomb gauge as this was done in many previous works |5J, |6[ |7j, §, |10fl. In the zero tem- 
perature limit the singlet free energy defined above coincides with the well known static potential. 
In fact the calculations on the static potential in Ref. [jll|] are based on Eq. (|ET|). At finite tem- 
perature F\(r,T) gives information about in-medium modification of inter-quark forces and color 
screening. The singlet free energy F\ (r, T) as well as the zero temperature static potential is defined 
up to additive constant C which depends on the lattice spacing. Since the temperature is varied by 
changing the lattice spacing to compare the free energy calculated at different lattice spacings we 
normalized it at the smallest distance to the following form for the T = static potential 

0.385 1.263 

V{r) = — + -j-r. (1.2) 

r Yq 

The above definition gives a very good parameterization of the lattice data for the zero temperature 
static potential calculated in Ref. [p. 

In Fig. [I] we show the singlet free energy calculated on 16 3 x 4 and 16 3 x 6 lattices together 
with the above parameterization of the T = static potential. As one can see from the figure, 
F\ (r, T) is temperature independent at very small distances and coincides with the zero temperature 
potential as expected. At large distances the singlet free energy approaches a constant value. This 
can be interpreted as string breaking at low temperature and color screening at high temperatures. 
Note that the distance where the free energy effectively flattens off is decreasing with increasing 
temperatures. This is another indication of color screening. Although the calculations have been 
done on quite coarse lattices we see that the results on F\ (r, T) show a fairly good scaling with the 
lattice spacing. This is shown in Fig. ^| where the singlet free energies calculated on 16 3 x 4 and 
16 3 x 6 lattices are compared at temperatures T ~ T c . 

2. Renormalized Polyakov Loop 

The expectation value of the Polyakov loop (L(r)) = (TrW(r)) is the order parameter for the 
deconfining transition in pure gauge theories. In full QCD dynamical quarks break the relevant 



2 



Singlet Free Energies 



Konstantin Petrov 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 



-F^r.T) [GeV] 




0.95T. 
0.99T C 
1 .00T. 
1 .021°. 
1 .05T„ 



0.2 



0.4 



0.6 0.8 
r[fm] 



1.2 




e e © e s> ©@ ® ©©©©gewegeft > 



'0.90T C 
0.97T. 
1.01T° 
1 .28T' 
1 .67T° 
2.12T, 
2.51 T' 



0.6 0.8 
r[fm] 



1.2 



Figure 1: The singlet free energy calculated on 16 x 6 (left) and on 16 x 4 lattice for m q = 0. lm s at different 
temperatures. The thick line shows the parameterization of the zero temperature potential discussed in the 
text. 
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Figure 2: The comparison of the singlet free energy for m q = 0.1m s calculated on 16 3 x 6 and 16 3 x 4 
lattices at T ~ T c . 



Z(3) symmetry explicitly and it is no longer the order parameter. Still it remains an interesting 
quantity to study the deconfinement transition as it shows a rapid increase in the crossover region 
|ii and can be used to determine the transition temperature [[T], [l3|]. The Polyakov loop 
defined above strongly depends on the lattice cutoff and has no meaningful continuum limit. On 
the other hand a correlator of Polyakov loops is a physical quantity and corresponds to the color 
averaged free energy up to a normalization constant. It satisfies the cluster decomposition 



exp(- J F(r,r) +C )= l -{L{r)V (0))| r _ 



l^(0)}| : 



(2.1) 



The normalization constant can be fixed from the color singlet free energy. Moreover, at large dis- 
tances the color singlet free energy and the color averaged free energy approach the same constant 
Foo(T). Therefore, following Ref. [Q], we define the renormalized Polyakov loop as 
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Figure 3: The renormalized Polyakov in the vicinity of the transition calculated on 16 x 6 lattices for 
different quark masses. Also shown in the figure is the corresponding result from 16 3 x 4 lattice. 



L ren (T) = exp(-^l). (2.2) 

Our numerical results for the renormalized Polyakov loop for different quark masses and two lattice 
spacings are summarized in Fig. |3[ One can see from the figure that L ren (T) shows an almost 
universal behavior as function of T /T c for all quark masses studied by us, including the case of three 
degenerate flavors. This suggests, that in the region of the small quark masses we studied, the flavor 
and quark mass dependence of the deconfinement transition can be almost entirely understood in 
terms of the flavor and quark mass dependence of the transition temperature T c . Note that the 
results obtained on 16 3 x 4 lattice are in remarkably good agreement with the results obtained on 
16 3 x 6 lattices, indicating again that the cutoff effects are small. 

It is interesting to compare our results for the renormalized Polyakov loop with the calcula- 
tions for three degenerate flavors performed with Asqtad action ^ as well as with the two flavor 
calculations with p4 action at larger quark masses JT0[]. The calculations with Asqtad action have 
been done on lattices 12 3 x 4 and 12 3 x 6. Critical temperatures are 194(15)MeV for N t = 6 and 
199(8)MeV for A 7 , = 4, m q = 0.2wj[|l]. This is consistent with out own estimates from the max- 
imum of susceptibility of the renormalized Polyakov loop. The two flavor calculations with p4 
action were performed on 16 3 x 4 lattice and for quark mass of about 1.54m s with m s being the 
physical strange quark mass. This comparison is shown in Fig. |] where the renormalized Polyakov 
loops is plotted in a wider temperature range. We see some discrepancy in the transition region be- 
tween the results of our study and earlier Asqtad calculations. For N t = 4 it can easily be explained 
by the large error for the determination of the critical temperature - we can pick a temperature in 
the allowed region which makes curves lie almost on each other. The N t = 6 case shows more 
dramatic descrepancy and here, apart from even larger ambiguity for the critical temperature, we 
have to point out that Asqtad simulation are done on a much smaller volume. 

Most noticable deviations from the results of previous calculations are seen at higher temper- 
atures. These deviations may come from the fact that the parameterization of the zero temperature 
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Figure 4: The comparison of the renormalized Polyakov loop calculated with p4 and Asqtad actions on 
Nt = 4 lattices (left) and Nt = 6 lattices (right). For iVj = 4 we also shown the results from 2 flavor p4 
calculations. 
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potential given by Eq. ( |1.2| ) may not be appropriate for the small lattice spacings corresponding 
to high temperatures. Also the parameterization of the non-perturbative beta functions used in the 
present analysis is based on the analysis of ro for gauge coupling j3 =6/g 2 in the range between 
3.3 and 3.4. It may not be accurate for j8 values corresponding to high temperatures. This also 
could yield a discrepancy in the value of the renormalized Polyakov loop in the high temperature 
region. 



3. Conclusions 

We have calculated the singlet free energy of a static quark anti-quark pair in full QCD in the 
region of the small quark masses. We have found that this quantity shows little cut-off dependence 
and can be calculated reliably on relatively coarse lattices. The renormalized Polyakov loop derived 
from the large distance limit of the singlet free energy has been also calculated and compared with 
earlier results. We have found that the flavor and quark mass dependence of the renormalized 
Polyakov loop can be absorbed almost entirely in the flavor and quark mass dependence of the 
transition temperature. 
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